m1=2500; m2=320; k1=80000; k2=500000; b1 = 350; b2 = 15020;
num2=[-(m1*b2) -(m1*k2) 0 0];
den2=[(m1*m2) (m1*(b1+b2))+(m2*b1) (m1*(k1+k2))+(m2*k1)+(b1*b2) (b1*k2)+(b2*k1) k1*k2];
'G(s)2'
printsys(0.1*num2,den2)
t=0:0.01:50;
y2= step(0.1*num2,den2,t);
plot(t,y2)
grid
title('open-loop response to 0.1 m step disturbance (transfer function)')